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I.  INTRODUCTION 
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Since  the  mid-1970s,  when  it  was  recognized  that  heavily  ionizing 
energetic  charged  particles  were  the  cause  of  logic  and  memory  upsets  in 
integrated  circuits,  these  "single-event  upsets"  have  been  the  subject  of 
considerable  research  and  analysis.  The  mechanism  for  producing  single-event 
upsets  is  well  understood.  Energetic  heavy  ions  can  cause  these  upsets 
because  they  can  deposit  more  energy,  via  ionization,  than  the  energy  of  the 
signal  representing  the  information  stored  in  a  logic  cell  of  an  integrated 
circuit.  Much  work  has  been  devoted  to  characterizing  the  sensitivity  of 
various  specific  integrated  circuits  and  logic  families,  and  to  determining 
methods  for  estimating  the  rate  at  which  a  given  device  may  be  upset  by  the 
radiation  environment  in  which  it  operates. 


Tests  to  determine  the  sensitivity  of  a  given  device  seek  to  measure  two 
characteristics:  (1)  a  threshold  energy  deposit,  or  critical  charge,  which  is 
the  minimum  amount  of  energy  that  must  be  deposited  in  a  sensitive  region  to 
upset  the  device;  and  (2)  the  size  of  the  sensitive  region.  (Deposited  energy 
and  charge  are  related  by  a  constant  of  proportionality  equal  to  3.6  eV  per 
charge  for  silicon.)  The  rate  at  which  a  device  may  be  upset  is  equal  to  the 
flux  of  particles  having  a  sufficient  ionization  and  pathlength  through  the 
sensitive  volume  such  that  the  product  (ionization  «  pathlength)  is  greater 
than  the  critical  charge.  Thus  the  device  upset  rate  depends  on  an  integral 
over  a  pathlength  distribution  of  the  spectrum  of  particle  ionization.  This 
particle  ionization  spectrum  is  also  called  the  linear  energy  transfer,  or 
LET,  spectrum.  The  LET  spectrum  is  determined  from  the  composition  and  energy 
spectra  of  all  species  of  ionizing  particles  in  the  environment,  together  with 
the  thickness  of  material  surrounding  a  given  device  which  shields  it  from 


this  radiation  environment. 
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Integrating  the  LET  spectrum  over  a  pathlength  distribution  is  a 
straightforward  task,  but  not  a  particularly  simple  one.  Recognizing  the  need 
for  a  "quick  and  dirty"  approach  to  this  problem,  Petersen  et  al .  (1983)  have 
proposed  a  "single-event  upset  figure  of  merit"  which  yields  an  estimated 
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upset  rate  for  any  device,  given  its  critical  charge  threshold  and  sensitive 
area.  This  "figure  of  merit"  calculation  incorporates  approximations  to  both 
the  geometry  of  the  sensitive  region  and  the  form  of  the  LET  spectrum,  so  that 
the  integral  of  (LET  x  pathlength)  takes  a  particularly  simple  form.  To 
derive  their  "figure-of -merit,"  Petersen  et  al.  adopted  an  approximation  to 
the  differential  flux  intensity  of  the  LET  spectrum  by  assuming  that  the 
differential  particle  flux  of  particles  with  a  given  LET  was  proportional  to 
the  -3  power  of  the  LET  (equivalent  to  the  integral  flux,  F  -  A  x  (LET)-2,  for 
some  constant  A).  The  result  of  these  approximations  is  that  the  "Petersen 
figure  of  merit"  upset  rate  for  any  device  may  be  written  as 

2 

R  -  K(S/L  )  (approximate) 
c 

where  R  is  the  f igure-of-merit  upset  rate;  S  is  the  surface  area  of  the 
sensitive  region;  Lc  is  the  device  LET  threshold;  and  K  is  a  constant 
proportional  to  A,  the  intensity  of  the  LET  spectrum. 

Petersen  et  al.  proposed  that  their  figure  of  merit  be  based  on  a 
specific  model  of  the  heavy-ion  environment,  the  "90*  worst  case"  environment 
defined  by  Adams  et  al .  (1981).  As  the  name  implies,  this  model  was  defined 
so  that  the  natural  fluxes  of  heavy  ions  from  the  galactic  cosmic  rays  and 
from  solar  and  interplanetary  acceleration  events  could  be  expected  to  exceed 
the  model  for  no  more  than  10$  of  any  large  number  of  randomly  selected,  short 
periods  of  time. 

Given  this  environment,  Petersen  et  al.  approximated  the  model 
differential  LET  spectrum  [with  LET  measured  in  (MeV/g/cm2) ] : 

f (differential  flux)  =  5.8  x  10®  (LET)-®  particles/ (cm2  day ) (MeV/g/cm2) 
which  is  equivalent  to  the  integral  LET  spectrum: 


F  (integral  flux)  =  2.9  *  10®  (LET)-2  particles/ (cm2  day) 


and  determined  that  the  constant  K  should  therefore  have  the  value 
K  ■  5  *  10*1®  (upsets/cell/day)  *  ( pC/ynO^/ym^ 

p 

(for  S  in  ym  and  for  Lc  in  picocoulombs  per  ym).  (In  another  system  of 

p 

units,  the  constant  in  the  equation  for  F,  above,  takes  the  value  0.1*45  (m  -s 
sr)-1,  with  LET  expressed  in  MeV/ym.  These  are  the  units  used  in  the  tables 
attached  to  this  report.) 

Petersen  et  al.  specifically  note  that  this  use  of  a  reference 
environment  and  the  errors  inherent  in  approximating  the  LET  spectrum  as  a 
power  law  over  a  wide  range  of  LET  do  not  permit  accurate  predictions  of 
absolute  upset  rates  to  be  made.  However,  the  results  of  their  figure-of- 
merit  calculation  may  be  useful  in  predicting  relative  upset  rates  of  various 
devices,  given  estimates  or  measurements  for  S  and  Lc. 

The  Petersen  figure  of  merit  is  appealing  because  of  its  simplicity. 

Thus  there  is  a  significant  impetus  to  apply  this  method  to  other  reference 
environments.  In  the  rest  of  this  report,  a  number  of  other  environment 
models  are  described.  LET  spectra  are  calculated  for  each  of  these 
environments,  and  revised  values  for  K  are  defined,  as  appropriate,  for  each 
environment  model. 


II.  REFERENCE  ENVIRONMENTS  FOR  THE  CALCULATION  OF  HEAVY-ION  LET  SPECTRA 


The  major  sources  of  energetic,  heavily  ionizing  particles  are  the 
galactic  cosmic  rays  and  energetic  ions  accelerated  during  solar  flare 
events.  Time-intensity  profiles,  energy  spectra,  and  the  composition  of  these 
two  populations  of  particles  are  sufficiently  different  to  warrant  separate 
treatment.  Galactic  cosmic  rays  are  continually  present;  their  intensities 
vary  relatively  slowly  over  the  11-year  solar  cycle,  declining  as  solar 
activity  increases  and  recovering  once  again  during  periods  of  minimum  solar 
activity.  In  contrast,  energetic  solar  ions  appear  near  earth  as  the  result 
of  impulsive  solar  flare  events.  Their  intensities  vary  rapidly  over  periods 
of  a  few  hours  to  days,  and  the  peak  intensity  of  the  heavy-ion  flux  from  a 
single  event  varies  greatly  from  flare  to  flare. 

Adams  et  al.  (1981)  have  made  a  comprehensive  summary  of  the  results  of 
several  decades  of  research  on  the  fluxes  of  the  galactic  cosmic  rays  and  of 
solar  flare  particle  events.  From  this  summary  they  have  derived  a  set  of 
model  environments  that  describe  the  interplanetary  heavy-ion  flux  and  fluence 
distributions  that  may  be  expected  under  a  wide  variety  of  conditions.  These 
models  have  since  become  standards  for  use  in  the  specification  of  system 
performance  and  immunity  from  single-event  upsets. 

The  data  base  of  energetic  heavy-ion  events  that  was  available  to  Adams 
et  al .  was  severely  limited,  however,  and  Chenette  and  Dietrich  (1984)  have 
published  a  summary  of  solar  flare  event  characteristics  over  the  solar  cycle 
from  1973  to  1984,  in  an  attempt  to  fill  this  gap. 

The  following  environment  models  have  been  selected  to  expand  the 
Petersen  figure  of  merit  method: 

1.  "Galactic":  the  solar  minimum  galactic  cosmic  ray  spectrum  described 
by  Adams  et  al .  ( 1 98 1 ) . 

2.  "IMP-8":  a  solar  flare  heavy-ion  spectrum  based  on  the  24  September, 
1977  solar  flare  event  as  described  by  Chenette  and  Dietrich  (1984), 
and  measured  by  an  experiment  aboard  the  satellite  IMP-8, 

3.  "Ordinary":  the  Adams  et  al .  "ordinary"  solar  flare  model, 


4.  the  "90J  worst  case"  flare  environment  model  of  Adams  et  al . ,  and 

5.  "AL":  the  anomalously  large  flare  model  of  Adams  et  al. 

Each  of  these  model  environments  is  described  in  terms  of  both  the 
instantaneous  peak  particle  flux  and  the  total  fluence  of  particles  expected 
(flux  integrated  over  time).  In  the  case  of  the  galactic  cosmic  ray  model, 
the  total  annual  fluence  during  a  year  near  solar  minimum  conditions  is 
compared  to  the  total  fluence  for  a  single  solar  flare  event.  Additionally, 
for  the  solar  flare  models,  two  kinds  of  particle  composition  are  studied: 
"normal"  composition  and  heavy-ion-enriched  composition  (see  the  descriptions 
of  each  model  for  details). 

Figures  1  through  4  and  Tables  1  through  4  summarize  the  integral  LET 
spectra  that  result  from  these  models.  In  each  case  the  LET  spectrum  was 
calculated  assuming  a  passive  shield  of  1  g/cm2  (equivalent  to  150  mils  of 
aluminum).  Figures  1  and  2  and  Tables  1  and  2  give  integral  LET  spectra 
corresponding  to  the  peak  event  flux  for  normal  (1)  and  heavy-ion-enriched  (2) 
composition.  Figures  3  and  4  and  Tables  3  and  4  give  integral  LET  fluence 
spectra  for  each  event,  assuming  normal  (3)  and  heavy-ion-enriched  (4) 
composition.  (Heavy-ion  enrichment  does  not  apply  to  the  galactic  cosmic  ray 
environment  models.)  The  results  at  LET  <  8  MeV/ym  are  poorer  approximations 
to  the  true  LET  spectrum  expected  in  space.  This  is  because  (1)  the 
abundances  of  the  elements  that  have  such  large  LETs  are  only  poorly  known, 
and  (2)  at  tnese  large  LETs,  as  a  result  of  the  small  fluxes  and  fluences, 
other  factors  (e.g.,  nuclear  interactions)  may  dominate  the  incident  primary- 
particle  flux. 

Additionally,  in  Tables  5  through  8  the  integral  LET  spectra  in  Tables  1 
through  4  are  presented  as  ratios  to  the  Petersen  environment.  The  values  in 
these  tables  may  be  used  to  multiply  the  Petersen  environment  directly  to 
obtain  the  LET  spectrum  of  each  of  Tables  1  through  4.  For  the  fluence  models 
reported  in  Tables  6  and  8,  the  multiplier  is  divided  by  1  day  (86,400  sec). 
The  results  that  appear  in  the  tables  may  then  be  read  as  an  equivalent  number 
of  days  of  exposure  to  the  Petersen  environment  or  as  a  multiplier  to  the 
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Peak  Event  Flux  LET  Spectra  for  Several  Solar  Flare 
Models  Assuming  Normal  Composition,  and  the  Galactic 
Cosmic  Ray  (GCR)  Model  LET  Spectrum  at  Solar  Minimum 
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Fig.  3.  Total  Event  Fluence  LET  Spectra  for  Several  Solar 

Flare  Models  Assuming  Normal  Solar  Composition,  and 
the  Annual  Galactic  Cosmic  Ray  (GCR)  LET  Fluence 
Spectrum  during  Solar  Minimum 
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Total  Event  Fluence  LET  Spectra  for  Several  Solar 
Flare  Models  Assuming  Heavy-Ion  Enrichment,  and  the 
Annual  Galactic  Cosmic  Ray  (GCR)  Model  LET  Fluence 
Spectrum  during  Solar  Minimum 
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Table  3.  Heavy- Ion  LET  Spectra  for  several  environment  models. 

Normal  composition  -  event  fluenee  -  1  g/cm^  shielding 
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daily  count  of  upset  rates,  calculated  from  the  Petersen  method,  that  may  be 
expected  during  the  course  of  the  event  (for  the  galactic  flux  an  "event"  is 
one  day  of  exposure). 
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III.  APPLICATION  OF  MODEL  LET  SPECTRA  TO  PETERSEN'S 
FIGURE-OF-MERIT  CALCULATION 


To  apply  the  integral  LET  spectra  of  the  previous  section  to  the  Petersen 
f igure-of-merit  calculation,  these  spectra  must  be  approximated  as 

F  -  A  *  (LET)-2 

Because  the  spectra  are  not  well  represented  by  such  a  form,  the  constant  A 
may  be  selected  in  a  number  of  different  ways.  If  the  LET  spectra  were  well 
represented  by  Petersen's  assumption  of  a  power-law  with  index  -  -2,  then  the 
entries  in  each  column  of  Tables  5  through  8  would  be  independent  of  energy. 
The  fact  that  they  are  not  illustrates  the  errors  inherent  in  the 
approximation. 

The  conservative  method  is  to  determine  the  value  of  A  so  that  the 
approximation  is  equal  to  or  greater  than  the  model  LET  spectrum  at  every 
point.  This  is  equivalent  to  choosing  the  maximum  value  in  each  column  of 
Tables  5  through  8,  and  this  is  the  choice  selected  in  this  report.  For  any 
choice  of  A  smaller  than  this  one,  it  is  possible  to  construct  an  example 
wherein  the  calculated  f igure-of-merit  upset  rate  would  be  smaller  than  the 
rate  calculated  from  the  proper  pathlength  distribution  integral.  Thus, 
without  ji  priori  knowledge  of  the  types  of  devices  to  which  this  method  will 
be  applied,  this  conservative  choice  is  the  best.  Table  9  presents  a  summary 
of  the  results  of  the  tables  and  compares  them  to  other  work. 


IV.  COMPARISON  TO  THE  RESULTS  REPORTED  BY  E.  C.  SMITH  OF  TRW 


In  a  private  report  of  limited  circulation,  E.  C.  Smith  from  TRW  has 
determined  a  set  of  multipliers  that  differ  significantly  from  those  derived 
in  the  present  report.  In  this  section  these  differences  are  explored  and 
explained. 

Smith  examined  four  cases,  all  of  which  assume  a  passive  shielding 

equivalent  to  1  g/cm2  of  aluminum:  (1)  a  solar  cycle  average  environment,  (2) 

an  "anomalously  large"  flare  without  any  heavy-ion  enrichment,  and  two 
"ordinary"-sized  flares  that  assume  (3)  a  "mean"  composition  and  (4)  a  " 90 % 
worst-case"  composition  (i.e.,  a  flare  with  a  significant  enrichment  of  heavy 
ions).  The  results  of  Smith’s  calculations  are  also  summarized  in  Table  9. 

One  of  the  consistent  differences  between  Smith's  results  and  the  results 
reported  here  is  that  Smith  chose  to  determine  the  constant  A  in  the  power-law 

approximation  to  the  LET  spectrum  by  using  a  weighted  average  method.  (In 

effect,  Smith  averaged  his  equivalents  of  the  entries  in  each  of  Tables  5 
through  8  over  some  energy  range.)  This  procedure  produces  constants  that  are 
smaller  than  the  upper-limit  conservative  choice  by  a  factor  ranging  from  1.5 
to  <3.  This  is  all  that  is  necessary  to  explain  the  differences  between 
Smith's  results  for  the  galactic  flux  models  and  those  of  Tables  1  through  8. 

Whether  or  not  this  averaging  technique  is  proper  or  appropriate  is  a 
matter  of  opinion.  It  is  my  opinion  that  it  is  not.  I  believe  that  when  such 
approximations  are  used  to  simplify  a  problem  like  this  one,  those 
approximations  should  yield  error  rate  estimates  that  are  upper  limits.  If 
the  upper  limit  error  rate  is  intolerable,  then  a  more  accurate  calculation 
can  be  made  to  refine  the  result.  If  the  approximate  result  is  not  based  on 
an  upper  limit,  then  the  actual  error  rate  may  be  larger  than  expected,  which 
may  lead  to  unknown  and  unaccountable  problems  in  the  effectiveness  of  the 
system. 

Another  correction  that  Smith  applied  to  his  results  for  the  galactic 
flux  environment  was  to  divide  the  solar  minimum  environment  model  result  by 


1.5  as  an  estimate  for  a  solar  cycle  average.  This  is  approximately  correct, 
but  it  is  important  to  remember  that  it  is  appropriate  only  for  solar  cycle 
(10  or  more  years)  averages. 

For  the  solar  flare  models  there  are  several  more  dramatic  differences 
between  Smith's  results  and  those  presented  here.  The  difference  of  a  factor 
of  1.5  between  Smith's  "AL  -  Mean  Composition"  multiplier  (3300.0)  and  the 
value  derived  here  (5000.0)  is  due  to  the  selection  effect  described  above. 
Just  as  in  the  galactic  cosmic-ray  flux  multipliers,  one  part  of  the 
explanation  for  the  differences  between  these  flare  multipliers  may  be  traced 
to  the  choice  of  how  to  fit  the  LET  spectra.  These  differences  lead  to  a 
consistent  bias  between  the  results.  There  are,  however,  additional  effects. 

When  comparing  results  for  the  flare  fluence  models,  the  multipliers  are 
dependent  not  only  on  the  spectrum  fitting  technique,  but  also  on  assumptions 
concerning  the  time-intensity  profile  of  the  flare.  For  the  IMP-8  flare  the 
total  fluence  was  measured  directly,  and  from  that  measurement  a  time  factor 
(equal  to  20  hr)  was  derived.  This  time  factor  is  the  ratio  of  the  total 
fluence  divided  by  the  peak  event  flux.  In  contrast,  the  Adams  models  provide 
separate  spectra  for  peak  flux  and  event  fluence  for  each  of  the  flare 
cases.  To  obtain  his  result  for  an  anamolously  large  flare  model,  Smith 
assumed  a  time  factor  of  11.08  hr.  The  justification  for  this  choice  is  not 
documented,  but  compared  to  the  Adams  model  result  it  produces  a  spectrum 
smaller  by  a  factor  of  2.6  (1540  days  rather  than  4000  days  equivalent 
exposure  relative  to  the  Petersen  environment). 

The  most  serious  problem  with  Smith's  analysis  concerns  the  corrections 
that  he  applies  to  his  AL  flare  multiplier  to  account  for  deficiencies  that  he 
finds  in  the  models.  In  an  appendix,  Smith  declares  that  the  model 
overestimates  the  proton  fluence  observed  in  the  August,  1972  flare  by  a 
factor  of  1.6  and  that  the  heavy-ion  abundance  was  overestimated  by  a  factor 
of  2.3.  Together  these  factors  serve  to  reduce  Smith's  result  from  1540  days 
to  420  days.  Combined  with  the  factor  of  2.7  difference  cited  above,  these 
corrections  lead  Smith  to  a  final  result  that  is  only  10$  of  the  Adams  model 
prediction. 
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While  Smith  may  argue  that  his  corrections  can  be  justified  by  a  specific 
event,  it  is  not  clear  that  they  are  appropriate  in  a  generic  model.  This 
conflict  raises  a  policy  question  like  that  of  fitting  the  LET  spectrum. 

Should  one  adopt  a  model  that  has  been  carefully  minimized  to  suit  a 
particular  function,  or  should  one  retain  a  standard  model  (like  that  of  Adams 
et  al.)?  It  is  important  to  remember  that  anomalously  large  flares  are 
extremely  rare,  and  that  the  data  base  that  documents  the  features  of  these 
events  is  extremely  limited.  The  uncertainties  of  any  parameters  that  3eek  to 
describe  AL  flares  as  a  class  are  large.  The  factor  of  10  reduction  in 
Smith's  value  for  the  "AL  flare"  multiplier  places  it  at  the  low  end  of  the 
range.  The  significant  concern  is  that,  by  adopting  an  optimistic  model  of 
the  environment,  real  problems  may  be  ignored  rather  than  solved. 


V.  RELATIONSHIP  OF  RESULTS  TO  OTHER  ENVIRONMENT  SPECIFICATIONS 


The  Adams  models  have  received  wide  circulation.  As  a  result,  they  have 
served  as  the  basis  for  a  number  of  working  specifications  for  model 
environments.  When  these  other  specifications  share  the  same  spectral  form  as 
the  Adams  models,  it  is  easy  to  convert  the  results  of  this  report  to  apply  to 
another  specification.  For  example,  if  one  built  a  model  assuming  a  10-year 
environment  including  five  AL  flares  -  one  at  the  intensity  of  the  August, 

1972  event,  two  at  one-half  that  intensity,  and  two  at  one-quarter  that 
intensity  -,  then  since  this  is  equivalent  to  a  total  of  2.5  times  the  August, 
1972  flare  fluence,  the  total  number  of  upsets  from  these  events  can  be 
obtained  from  a  fluence  spectrum  that  is  2.5  times  greater  that  that  of  Tables 
3  or  il.  Similarly,  other  combinations  of  these  standard  galactic  flux  and 
flare  models  can  be  constructed. 


